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The number of conducting channels for one propagating direction is equal to that for the 
other direction in ordinary quantum wires. However, they can be imbalanced in graphene 
nanoribbons with zigzag edges. Employing the model system in which a degree of channel- 
number imbalance can be controlled, we calculate the shot-noise power at zero frequency 
by using the Boltzmann-Langevin approach. The shot-noise power in an ordinary diffusive 
conductor is one-third of the Poisson value. We show that with increasing the degree of 
channel- number imbalance, the universal one-third suppression breaks down and a highly 
nonuniversal behavior of shot noise appears. 

KEYWORDS: perfectly conducting channel, graphene nanoribbon, Boltzmann-Langevin equa- 
tion, one-third suppression 

1. Introduction 

More than a decade ago, Barnes et al. 1 ' 2 ' considered a fictitious disordered quantum wire 
of length L having N conducting channels in minority direction and N+m conducting channels 
in majority direction (m = 0, 1, 2, ... ), and showed by using a scattering argument that this 
system has m perfectly conducting channels in the majority direction. This indicates that if 
m > 1, the conductance in the minority direction vanishes in the limit of L — > oo while that 
in the majority direction converges to (e 2 /27r)m. That is, Anderson localization arises only in 
the minority direction. Although such an intriguing behavior was predicted, its detail has not 
been investigated until recently. This is probably because a realistic example of the system 
with channel-number imbalance has been lacking. 

A few years ago, Wakabayashi et al. 3 ~^ have pointed out that graphene nanoribbons with 
zigzag edges are a promising candidate for the system with channel-number imbalance. The 
band structure of a zigzag nanoribbon has two energy valleys, called K_ and K + valleys, well 
separated in momentum space, and each valley has an excess one-way channel arising from 
a partially flat band. 6 -* Due to the presence of an excess one-way channel, the numbers of 
conducting channels for two propagating directions are imbalanced in each valley. If impurity 
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Fig. 1. The band structure of a quantum wire with two energy valleys K- and K + , in each of which 
the numbers of conducting channels are imbalanced between two propagating directions. 



potentials are long-ranged and thus intervalley scattering is absent, the two energy valleys are 
disconnected leading to the realization of a channel-number imbalanced system. Inspired by 
this observation, the statistical behavior of the conductance in disordered wires with channel- 
number imbalance has been studied extensively. 7 ~ 14 ) Recently, the present author studied the 
averaged conductance of this system in the incoherent regime, where the phase coherence of 
conduction electrons are completely lost, and showed that its characteristic behavior reflecting 
the presence of perfectly conducting channels is unchanged even in the absence of phase 
coherence. 15 ^ 

In this paper we study nonequilibrium shot noise in the system with channel-number 
imbalance. In an ordinary diffusive conductor with conductance G, the shot-noise power P s hot 
at a bias voltage V is given by P s hot = (2/3)e|y|G, 16 ' 17 ) which is one-third of the Poisson 
value -Ppoisson = 2e|V|G. This behavior, known as the one-third suppression of shot noise, has 
been confirmed in several experiments. 18 ~ 20 ) The one-third suppression arises in the diffusive 
regime regardless of the presence or absence of the phase coherence as long as quantum wires 
are shorter than energy relaxation length. 21 ) Does this universality hold even in the presence 
of channel-number imbalance? To answer this question we employ the model system in which 
a degree of channel-number imbalance can be controlled, and calculate -P s hot m the incoherent 
regime. We show that with increasing the degree of channel-number imbalance, the one-third 
suppression breaks down and a highly nonuniversal behavior of shot noise appears. We set 
h = &b = 1 and ignore the spin degeneracy throughout this paper. 

2. Model and Formulation 

We describe our model system for disordered quantum wires with channel-number im- 
balance. 15 ) As schematically shown in Fig. 1, its band structure has two energy valleys, K- 
and K + . We assume that the numbers of conducting channels are imbalanced in each valley 
between two propagating directions although the total number of conducting channels in one 
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direction is equal to that in the other direction. For concreteness we assume that the K_ val- 
ley has N + to right-moving channels and N left- moving channels, and the K + valley has N 
right-moving channels and N + m left- moving channels (see Table I). Figure 1 corresponds to 



Table I. The number jVr of right-moving channels and the number 7Vl of left-moving channels in 
each valley. 



Valley 


N R 


N L 


K- 


N + m 


N 


K+ 


N 


N + TO 



the case of N = 2 and m = 1. When m > 1, a degree of channel- number imbalance can be con- 
trolled by the strength of intervalley scattering due to impurities. Apparent channel-number 
imbalance appears if intervalley scattering is absent, while it disappears with increasing the 
strength of intervalley scattering. Let ( e Lnfc) ^ e ^ ne ener gy °f an electron state with wave 
number k in the nth right-moving (left-moving) channel of the K± valley. We introduce the 
corresponding distribution function g z Xnk -, where X = R, L and z = +,—. The group velocity 
is given by v z Xnk = de z Xnk /dk. 

We assume that our model system of length L is placed along the x axis, and is connected to 
left and right reservoirs at x = —L/2 and L/2, respectively. A bias voltage V is applied between 
the two reservoirs, and thus the electron distribution is described by /l(e) = /fd(e — e ^/2) 
in the left reservoir and /r(£) = /fd(e + eV/2) in the right reservoir, where /fd represents 
the Fermi-Dirac function. We study noise power of this model on the basis of the Boltzmann- 
Langevin approach. An expression of noise power has been derived in refs. 22 and 23 for 
disordered conductors in a continuum representation. To adapt it to our model we must take 
account of its unique band structure by introducing the concept of subbands. We outline the 
derivation of noise power in a subband representation. 

We introduce rj = {Xnkz} to characterize electron states. For later convenience we rep- 
resent the time-reversed counterpart of rj as fj = {Xnkz}, where R = L, L = R, k = —k, and 
i = =F. We decompose g v (x,t) into the average and fluctuations as 

g v (x,t) = f v (x) + Sf v (x,t), (1) 

where f v (x) = {g r ,{x, t)) with (• • • ) indicating the average over t. The averaged part obeys the 
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Boltzmann equation in a subband representation, 24 ) 

( V r,-^ ~ limp®) fr,(x) = 0, (2) 

under the boundary condition f v (—L/2) = and f v (L/2) = fR(e v ). The collision integral 

for impurity scattering is given by 

/imp ® f V (x) = W v rf (f v > (x) - f v (x)) (3) 

n> 

with Wqtf being the scattering probability between the state with rf = {X'n'k'z 1 } and that 
with i] = {Xnkz}. The scattering probability is expressed as 

W v rf = 2-kM^S (e v > - e v ) (4) 

with Mjjrf = (| Ujj^i | 2 )imp; where U v ^i is the matrix element of impurity potential and (• • • )j mp 
indicates an ensemble average. The time-dependent part obeys 

(§i + Vv ~iL ~ /imp ®) 6 fv( x > *) = $M x i *)> ( 5 ) 
where 8j n represents the source of stochastic fluctuations. This has zero average, (Sj r) (x, t)) = 
0, and satisfies 25 ^ 

<c^(x,t)tyV(x / ,t')> = L5(x - x')5(t - (6) 

with 

J V:V '{x) = $ v , v >^2w V:V/ > [f v (l - f rj ») + f v »(l - f v )] 

rf 

-/„)], (7) 

where f v = f n (x) . The average and fluctuations of the current are expressed as 

(I)(x) = jJ2 v M^ (8) 
v 

5I(x,t) = ^^2v v 5f v (x,t), (9) 

v 

respectively. 

It is convenient to introduce the Green's function Q V;r) /(x,x';t) which satisfies 
( d d \ , , 

( ^ + Vr i-fa ~ J im P ® J Gri, v '(x,x';t) = 8 r) ^8{x - x')5(t). (10) 
We can express 5f v as 

8f v (x,t) = f dt' I dx'YsGwiXix'-t-t'^ix'j) 

J —CO J i 
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E f dt ' E G^(x,x'-t-t')v v ,5f vl {a(^ + 5),t'), (11) 



er=± ' 



7j'((7 V <0) 



where 5 is a positive infinitesimal. The second term is added to represent incoming fluctuations 
from the left and right reservoirs. In the reservoirs near the sample region, the incoming 
fluctuations obey 

(5f v (x,t)5f v ,(x',t')} = L5{x -x'- v v (t - 0)<W/x(l " fx), (12) 
where fx = /l^) for x,x' < —L/2 and fx = /r^) fo r x,x' > L/2. We define T v as 

roc 

T v (x) = dt^2v v >g- lfj (xo,x;t) + 6(x -x), (13) 

r;' 

where 9{x) is the Heaviside step function and we set xo = —L/2. This represents the trans- 
mission probability that an electron at x with fj arrives at xq = —L/2 and then goes into the 
left reservoir. Alternatively, noting the presence of time-reversal symmetry, we can interpret 
it as the transmission probability that an electron injected into the quantum wire from the 
left reservoir arrives at x with i]. We can show that T v {x) satisfies eq. (2). The boundary 
condition is given by 

Ti nk {-L/2) = 1, (14) 

T£ nk (L/2) = 0. (15) 

In terms of T v (x), the averaged part of the distribution is expressed as 

f v (x) = T„(s)ft( e „) + (1 - T v (x))f R (e v ). (16) 

By using eq. (16), the conductance defined by G = limv^o{I) (x) /V is obtained as 

G = e ^Y. v Mx)(-^(e r ,)Y (17) 
^ \ / 

which is independent of x due to the continuity of current. We consider the noise power at 

zero frequency 

/oo 
dt{5I(x,t)5I(x,0)). (18) 
-oo 

This is also independent of x, so we set x = xq without loss of generality. Using eqs. (6)-(16), 
we can express P in terms of T v (x). After manipulation (see ref. 23 for details), we obtain 

P = ^shot + Pth With 

Pshot = \\ d^^^^W " T-,{x))\{x){l - T v ,(x))(f L - f R ) 2 , (19) 
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^h = ^£ vr, 2T v (~)f L (l - f L ) - T ri (-^f (l - [/ L (l - /l) - /r(1 - /r) 

(20) 

where /x = fxi^ri)- Here -P s hot and P t h represent the shot-noise power and the thermal noise 
power, respectively. 

3. Calculation of Noise Power 

We hereafter ignore the /c-dependence of T Xnk (x) because transport properties are deter- 
mined by its value at the Fermi level. For simplicity we assume that v^ nk = and v £ nfc = 
with v = = — > 0. Furthermore, we assume that the scattering probability is determined 
by only the valley indexes z and z', and does not depend on details of initial and final states. 
Hence M z ^ k X , nlk , = M z > z ' . We define k z ' z ' = LM Z > Z ' /v 2 with k ±j± = k and = k' . Note 
that k and k' describe intravalley scattering and intervalley scattering, respectively. The trans- 
mission probability T Xn {x) does not depend on n due to the assumptions presented above. In 
terms of the mean free path lx( x ) presented in ref. 15, we can express 

T Xn (x)= 1 --^ + si g n(v x ) l ^l. (21) 

The mean free path is given as l+(x) = l^(x) = l^(—x) and l-(x) = l^(x) = l^(—x) with 

=F m(K — k') (c(x) — cq) 



2 



l±(x) = x +2 + z 



m i( K 2 _ K ' 2 ) l 

+ K —= (d(x) - d ) - 2(2N + m)qK'c (x + -) 

V a 2 . 

where a = (k + K')[m 2 K + (8iV 2 + 8Nm + m 2 )/«'], q = (2N + m)(K + k'), and 



(22) 



c(x) = —t= cosh(y/ax) — sinh(y / ax), (23) 
V a 

d(x) = -= smh(y/ax) + cosh(y / ax). (24) 

The remaining constants are given as cq = c(— L/2), do = d(—L/2), and 

t 2/2 

S = 4(2iV + mW(l + %-] c + 2m 2 K ~^ d . (25) 
V 2 / y'a 

We can show that eq. (21) satisfies eq. (2) under the boundary condition of eqs. (14) and (15). 
The difference between and Z_(x) should be regarded as a manifestation of channel- 

number imbalance. 

Substituting eq. (21) into eq. (17) we obtain 

e 2 2ac 

° ~ 27 ( K + K ')£ ( } 

As shown in ref. 15, the averaged conductance G rapidly converges toward (e 2 /27r)m with 
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increasing L at k' /k = 0. If k'/k is sufficiently small but finite, G rapidly decreases to (e 2 /27r)m 
in an initial stage and then crosses over to a slow decrease toward zero. The slow decrease is 
accelerated with increasing k'/k. From eq. (20) we obtain = 4GT which is equivalent to the 
result for ordinary diffusive conductors. We next consider the shot-noise power. Equation (19) 
contains the summations over 77 = {Xnkz} and rf = {X'n'k'z'}. Carrying out the integrations 
over k and k' , we obtain 
e 2 Fv,T 



shot 



I dx E E ^ ( T l» - Ti nl {x)) 2 T* Xn (x) (l - z£ B ,(s)) , (27) 

J XnzX'n'z' 

where F V>T = e\V\ coth(e|V|/2T) - 2T. Note that Fy jT = e\V\ when e\V\ > T. We define 
the shot-noise factor p s hot = -Pshot /-Ppoisson with -Pp isson = IFy^G. This characterizes the 
suppression of shot noise and is equal to 1/3 in ordinary diffusive conductors. Substituting 
eq. (21) into eq. (27) and using eq. (26) we obtain 



(k + k')Y> 



2k{N + m)N{l + {x) + L(-x)) 2 L 2 + 4(i+(-x) + x) (/_ (x) — x) 

+ K'(iV + m) 2 (/+(x) + /+(-x)) 2 [L 2 + 4(/ + (x) - x) (l+(-x) + x) 
+ k'N 2 (/_ (x) + Z_ (-x)) 2 L 2 + 4(Z_ (x) - x) x) + x) 



+ 2k' (N + m)iV(; + (x) - L(x)) 2 I? - 4(Z+( -x) + x) (/_(— x) + x) 



• (28) 



We briefly examine two limiting cases without channel-number imbalance. One is m = 
with arbitrary k'/k and the other is k'/k = 1 with arbitrary m. The mean free path is expressed 
as l±(x) = (x+L/2)/A with A = 1 + N(k+k')L in the first case and A = l + (2N+m)KL in the 
second case. It should be noted that l+(x) = l-(x) indicates the absence of channel- number 
imbalance. In both the cases, we obtain 



Pshot = - (1 - A 3 ) 



(29) 



when A > 1, indicating the universal one-third suppression of shot noise. We present the 
result for p s hot in the presence of channel-channel imbalance for the case of N = 10. Figure 2 
displays p s hot a s a function of kL in the no inter- valley scattering limit of k'/k = for m = 0, 
1, and 2. We observe that p s hot is strongly suppressed below the universal value 1/3 and the 
suppression becomes large with increasing m. This is attributed to the fact that G rapidly 
converges toward (e 2 /27r)m with increasing kL while P s hot monotonically decreases. Figure 3 
displays p s hot as a function of (k + k')L at m = 1 for k'/k = 10°, 10~ 2 10~ 3 , 10 -4 , and 
10 -6 . The behavior of p s hot is very anomalous when k'/k < 1. As seen at k'/k = 10~ 3 , after 
an initial increase in the quasi-ballistic regime, p s hot decreases with increasing (k + k')L as 
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Fig. 3. p sho t as a function of (k + k')L when TV = 10 and m = 1 for k'/k = 10°, 10~ 2 , 10~ 3 , 10" 4 , 
and 10- 6 . 



in the case of k'/k = 0, but then turns to increase toward the universal value of 1/3. This 
behavior reflects the fact that channel-number imbalance becomes irrelevant in the long-L 
limit unless k'/k = 0. Indeed, the mean free path in this limit is approximately given by 
l + = l_ = (x + L/2)/A with 

A-1 + (2JV + '" ) " K ' J ,, L (30) 

2{2N + m)K' + m 2 ^- 

except for the vicinity of x = ±L/2 and hence the influence of channel- number imbalance 
disappears. This is consistent with the observation given below eq. (26). At k'/k = 0, the 
difference between l + and I- never vanish and therefore channel-number imbalance remains 
to be essential even in the long-L limit. 
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4. Summary 

We have calculated the shot-noise factor p s hot j which characterizes the suppression of shot 
noise compared with the Poisson value, in disordered quantum wires with channel-number 
imbalance in the incoherent regime. We have employed the simple model in which a degree 
of channel-number imbalance can be controlled by the relative strength k'/k of intervalley 
scattering. Channel-number imbalance completely disappears at k'/k = 1, while its influence 
becomes significant with decreasing k'/k. We have shown for m > 1 that p s hot decreases with 
increasing wire length L after an initial increase in the quasi-ballistic regime, but then turns to 
increase toward the universal value of 1/3 unless k'/k = 0. The suppression of p s hot becomes 
large with decreasing k'/k. This anomalous behavior should be contrasted to the universal 
one-third suppression of shot noise widely observed in ordinary diffusive conductors. 
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